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Abstract 



We are interested at first in the study of the monotonicity for the 
period function of the conservative equation (1) x+g(x) = 0. Some 
refinements of known criteria are brought. Moreover, we give neces- 
sary and sufficient conditions so that the analytic potential of equation 
(1) is isochronous. These conditions which are different from those in- 
troduced firstly by Koukles and Piskounov and thereafter by Urabe 
appear sometime to be easier to use. We then apply these results 
to produce families of isochronous potentials depending on many pa- 
rameters, some of them are news. Moreover, analytic isochronicity 
requirements of parametrized potentials will also be considered 
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1 Introduction and statement of results 

Consider the scalar equation with a center at the origin 




(1) 



or its planar equivalent system 




(2) 
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where x — x — and g(x) is analytic on R. 
Let be the potential of equation (1) 

G(x)= f "g(£)dt. 
Jo 

We suppose in the sequel the following hypothesis holds 



/_,s J There exist a < < b such that G(a) = G(b) = c, G(x) < c 
\ and xg(x) > for all a < x < b and i^O. 

Moreover, without lose of generality we will assume in the sequel that 
0(0) = and g'(0) = 1. 

Under these conditions system (2) admits a periodic orbit in the phase 
plane with energy c and a, b are the turning points of this orbit satisfying 
inequalities a < a < 0, < b < b, for some real a and b. 
This means the origin is a center of (2). This center is isochronous when 
the period of all orbits near G R 2 are constant (T = , 2?r = 2ir). The 

corresponding potential G(x) is also called isochronous. 

Let T(c) denotes the minimal period of this periodic orbit depending on 
the energy. It can be expressed under the following form 

T(c) = V2 f dX (3) 
Ja ^c-G{x) 



The period function is well defined for any c such that < c < c and when 
c < +oo then c is the energy of the homoclinic orbit with a < < b its 
turning points. It is well known that T(c) is an analytic function of c. 

Since the potential G(x) has a local minimum at 0, then we may consider 
an involution A by 

G(A(x)) = G(x) and A(x)x < 

for all x G [a, b]. So, any closed orbit is A-invariant and A exchanges the 
turning points: b = A(a). 

In fact, A(x) is well defined in the interval [a, b]. To see that, set the function 

x - A(x) 
P( x ) = o • 
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This function is such that p(A(x)) = —p(x) and p'(x) = — ■ Since 
A'(x) < we get p'(x) > and therefore p is an analytic diffeomorphism on 
[a, b}. Then A(x) = p^ 1 (—p(x)) is well defined. 

Conversely, by using this involution we may calculate for a prescribed pe- 
riod function T(c) the distance between the turning points. Indeed, Following 
Landau and Lifchitz we have [9, Chap. 3, 12.1]): 

which implies 

A(a)-a = b-a=^ r f ^2£. (4) 

When the potential (or the center 0) is isochronous or equivalently the 
period is constant for all orbits near zero one has necessarily 

b - a = 2V2c. 

The behavior of the period function is important in applied mathemat- 
ics. Many authors ([3], [4], [10], [11]) are interested in the monotonicity of this 
period function as well as in isochronicity cases of equation (1) ([6], [7], [10]). 
Isochronous centers plays a central role in several problems of dynamical 
systems. Also, it intervenes in the investigation of the quantum spectrum. 
Many studies have been devoted to the problem of the relationship between 
isochronicity equivalence and quantum isospectrality condition. See [1],[6] 
for more details. 

In many papers one tried to characterize isochronous potentials. In the case 
of the rational potentials, it can be shown [2] that the only rational potential 
with a constant period which not reduces to a polynomial is the isotonic one 

G{x) i n.r • 1 ' ' 



cr ax + 1 

The paper is organized as follows : 

- In the first part we study the monotonicity of the period function of a 
closed orbit of equation (1) depending on the energy and we produce better 
criteria than those known. In particular, we complete the work started by 
Chow and Wang [5]. 
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- We propose to analyze the notion of isochronicity in using a new track. 
We will provide necessary and sufficient condition so that the center of sys- 
tem (2) be isochronous, or equivalently its analytic potential is isochronous. 
Many applications are brought. 

- In the analytic case, let us write the potential under the following form 

G{x) = l -x 2 + G^x) + G 2 (x) 

where Gi(x) is odd and G 2 (x) is even. It is wellknown that the only symmet- 
ric isochronous potential is the harmonic one (i.e. G(x) = \x 2 ). We prove 
here that when G(x) is isochronous G\ and G 2 are related and consequently 
G 2 (x) = implies necessarily G\(x) = 0. 

Moreover, a simple proof of Chalykh-Veselov result [2] concerning rationals 
potentials will be provided. 

- Finally, we derive a three-parameters family of isochronous potential 
more general than the one given by Dorignac [6]. This family includes har- 
monic and isotonic and Urabe potentials as well as the Bolotin-Mac Kay 
potential. Other families of isochronous potentials also will be considered. 



2 The period function 

Let us consider the period function of equation (1) depending on the energy 
T = T(c). When g is an analytic function defined for x G [a,b], T(c) is 
also an analytic function defined for c G]0,c] such that lim c ^, +T(c) = 2n. 
The following result states sufficient conditions for T(c) to be monotonic. 

Theorem A Let g{x) be an analytic function and G(x) = f$ g(£,)d£, be 
the potential of equation (1) 'x + g(x) = 0. Suppose hypothesis (H) holds, 
let us define the n-polynomial with respect to G, 

f n {G) = /(0) + f(0)G + \f (0)G 2 + .... + ±f n (0)G n 

where /(0) = -g"(0)/3 = -(1/3)0(0), /'(0) = -(7/%'' 3 (0)+<?( 4 )(0)/5, 

/"(0) = -^ {6) (o) - fV 5 (0) + 2/ 2 (oy 4 )(o), ... 
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Suppose that for a fixed n <E N and for x G [0,6] one has 



(O 



> /»(<?) > ^[^(^))] 




Proof of Theorem A 



Recall that 



T(c) = V2 f 

J a 



■b 



dx 



In order to study the monotonicity of the period function T(c) depending on 
the energy, it is sometime convenient to study its derivatives. We need the 
following 

Lemma 2-1 The derivative of the period function (depending on the 
energy) T'(c) = ^ may be written 



This lemma has been initially proved by Chow and Wang [5] but the proof 
we give below is different. 

Proof Consider the change G(x) = s 2 , where s = u(x) is a function 
of x (this change has a sense because G(x) is positive in a neighborhood of 
according to the hypothesis xg(x) > 0). It yields 



By another change s = y/csin.9 the period function may be expressed 






T(c) = y/2f 



V2 / (u- 1 )'(y/cwie)dB 



72 u'(u^ 1 ( y /c) sin#)) 



J-n/2 
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Then we get by another way the derivative of the period function 
T'(c) = 4= {u~ 1 )"(Vcsme)sm6de. 

y/C J-tt/2 

By splitting T'(c) = X° 7r / 2 + "Jj Jo^ an< ^ usm g the formula 
(w -1 )"(w(s)) = — ("/(g))3 and making another time the change of variable 
s = u(x) = y/csin.9 we then obtain 



1 fb g 2 ( x ) -2G(x)g'(x) g(x 



dx. 



cyfiJo g\x) _ 

Since < the derivative can be written 

i rVW-aoWW^ 



To prove Theorem A observe at first 



T'(c) 



, 7^) / / =, = 9(x)dx. (6) 
cV2 •/« Jc-G(x^ 



dx. 



Write 

g 2 (x)Jc- G(x) 

By using the involution A and since x(g(x) > for rr 7^ we get 

cV2 ; A y/c-G(x) y cV2 J Jo ^J c -G(x) 

~ ( cv / 2 ) A y^T^) ^ + W io fiZc(xj 9{x)dX 
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- [ cV2 } Jo 9{X}dX - 

Then 

T (c) > (or <) — = / 2 == g{x)dx. 

cV2 J a ^c-G{x) 

Let us integrate by parts the following for p > 1 

g(x)dx = / dG = -G p V^G+2p / G^V^GdG. 

Jr-Cl(>A J Jr-rJ-A J 



c-G(x) J ^Jc-G(x) 

Therefore, it is easy to see that the following integral should be zero since 
G(a) = G(b) = c 

f b G p 

/ =g(x)dx = 



x 



for any integer < p < n. 

Thus, the right side of (6) is non zero and condition 

^-[G(x)/g 2 (x)]-f n (G)^0 
implies T'(c) ^ 0. This means T = T(c) is monotonic. 



As applications of Theorem A, we derive at first monotonicity conditions 
for the period function T(c) depending on the energy. This problem holds 
importance in dynamical systems (see [4], [10] and [11] for example). In this 
part, we find again some known criteria for the monotonicity of the period 
function and propose some others which are better and seem to be new in 
the literature. 

For a complete study and a comparison between these sufficients condi- 
tions we refeer to [3] and [4] and references therein. 

Notice that the monotonicity criteria for the period function produced by 
Cor. 2.5 of [5] is more general than those given by F. Rothe [10] and R. 
Schaaf [11]. 

We propose the following which slightly improves Cor. 2.5 of [5]. 
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Corollary 2-2 Suppose hypothesis (H) holds and let g(x) be an analytic 
function for a < x <b and 

G(x) = Jo g{£,)d£, be the potential of (1) and g"(0) ^ 0. 
1 - Suppose condition (C ) holds, this means 



9 2 {x) + ^j^g 3 {x) - 2G(x)g'(x) > 0(or < 0) 

then the period function of (1) is such that T'(c) > 0(or < 0) for < c < c. 
2 - Suppose 

±[G(x)/g 2 (x)} = a (7) 

where a = —g"(0)/3 ^ then is an isochronous center of (1). 
Moreover, the so-called Urabe potential 



, 4 2 „vT- ax. 
G i = - - (a + 2- 

where \ x \< 1/a is the unique analytic solution of (7) verifying 
g (Q) = G(0) = and G"(0) = 1. 



Proof of Corollary 2-1 Indeed, (7) may be written 

g\x) + 9 -^g\x) - 2G(x)g'(x) = g 3 (x)-^([G(x)/g 2 (x)] + 9 -^-g 3 (x). 
Then (7) must be 

d 



Implying 



(Jx [G(x)/<?(x)] = -gr(0)/3 = a 



G ( x ) = t;9 2 ( x ) + axg 2 (x) 



since g'(0) = 1. Thus, starting from (6) the derivative of the period function 
T(c) may be expressed under the following forms 

1 N [* g\x) - 2G(x)g'(x)_ dx 



T'( c ) = ( 1 ) f r{x)-2Ujx)g\, 

cV2 Ja gl( X )J C -G{x) 
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• 1 x f b U G l9 2 ) ( u 
— -=) I , =q{x)dx 

cV2 J Ja j c -G{x) 



= ( — = ) / , =g(x)dx. 
cV2 J Ja ^c-G(x) 

Integrating by parts, one gets 

r(c) = (-±=) f a = g(x)dx = 0. 
cV2 Ja ^Jc-G(x) 

On the other hand, since 

it implies that G = G(x) may be inverted and that x can be expressed 

x = y/2G + aG. (7') 

Remarks 2-3 Inversing now equation (7') and according to initial con- 
ditions g(0) = G(0) = and G"(0) = 1 it is easy to see that the potential 



„, . 4 2 a/1 - 
G(x) = — --(x + 2- 

is the unique solution of (7') or (7). G(x) is called Urabe potential. 
Notice that its derivative g(x) must also be a solution of 

5(g"(x)) 2 - Sg'(x)g"'(x) = 

and consequently for any Urabe potential we get necessarily 

S (4) (0) = f ^(0). (8) 

When in addition g"(0) = (i.e. a = 0) equation (7) has as unique 
solution the harmonic potential G(x) = (l/2):r 2 . 

In fact, as we will see below there exist isochronous potentials that do not 
verify the restrictive condition (8). Consider the following expression 

Hx) = ^(x) + ^ 3 (x)-2G(x)^) + (^^-^)^(x)G(x). (9) 



The following brings a new monotonicity condition for the period function 



Corollary 2-4 Let g(x) be an analytic function and 
G(x) = Jo g(£)d£ be the potential of equation (1). Suppose (C-J holds, this 
means 

9 2 {x) + 9 -^g\x) - 2G{x)g'{x) + _ ^M )g s {x)G{x) > Q(or < Q) 

for a < x <b then T'(c) > 0(or < 0) for < c < c. So a sufficient condition 
for (1) to have an isochronous center is 

Hx) = 9 2 {x) + 9 -^g\x) - 2G{x)g\x) + (J^- - 9 -^)g\x)G{x) = 
or equivalently 

±[G{x)/g\x)]=a + PG 



Proof of Corollary 2-4 Indeed, recall that 

^ = >^ = 

Then T{x) = means that 

9 2 - 2Gg' 



g 3 



= a + (3G 



where a = -g"(0)/S and f3 = -(7/9)a 3 + 7 /5 = -(7/%" 3 (0) + # (4) (0)/5. 
Consider again the derivative of the period function 



T'(c) = (— =) / f '=L g(x)dx = (—=) / y 'M= g{x)dx. 



So one obtains 

T'(c) = (^) f P==M x)dx 
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1 , f J a , . 1 , f fc (3G 

--) / , =g[x)dx + { — -=) I =g{x)ax. 



cV2 J* - G(x) ' cV2 J* yfc - G(x) 

By Corollary 2-2, the first integral of the right side should be 0. The second 
integral may be written as 



( — 7=) / i =g(x)dx = ( — 7=) / , =g(x)dx 

cV2 J Ja J c -G(x) y cV2>Ja ^J c _ G{x) 

= (-L) f ^ c - G{x)g{x)dx = (-L) / G(fc) /^v^GrfG = 

C\/2 Ja C\> 2 JG{a) 

Thus, one gets 

T'(c) = 0. 
Moreover, we also may deduce 

2 -^(^/9) = P9^G + ag/VG. 

On the other hands G = G(x) may be inverted and x can be expressed 
in terms of G. By integration and according to hypotheses G(0) = g(0) = 
and g'(0) = 1, it yields 

2(VG/g) = {2(3/3)GVG + 2aVG + V2, 

or equivalent ly 

l = (P/3)gG + ag + (-j^) 



So by another integration one obtains 

x= {f3 /6)G 2 + aG + V2G 

3 v 54 30 ; 

Remarks 2-5 Notice that clearly Corollary 2-4 is more general than 
Corollary 2-2. We then obtain a better criteria of the monotonicity for the 
period function. In the sense that (C ) implies (CJ which implies T'(c) > 
(or < 0). 

By the same way for any fixed n many other sufficient conditions ensuring 
the monotonicity of the period function of the form 

g 2 (x) - 2G(x)g'(x) - ag 3 (x) - (3g 3 (x)G(x) - -fg 3 (x)G 2 (x) - > (or < 0) 
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may be deduced. In particular, condition 

g 2 {x) - 2G{x)g'{x) - ag 3 {x) - (3g 3 {x)G{x) - ig 3 {x)G 2 {x) - = 

implies that the potential G(x) is isochronous and x can be expressed 

x = V2G + aG + (f3/6)G 2 + {i/24)G 3 + .... 

Thus, for a n-polynomial f n (G) = f {0) + f'(0)G + \f"(V)G 2 + .... + ±f n {0)G n 
conditions of Theorem A 

K) > «>/.(0 > 

ensuring T'(c) > (or < 0) are such that (C n+1 ) is finest than (C n ). 

3 Isochronicity conditions for a center of equa- 
tion (1) 

3.1 Isochronicity conditions for a center of equation (1) 

The problem to determine whether the center is isochronous has attracted 
many researchers for long time. This problem has been recently revived due 
to advancement of computer algebra. New powerfull algorithms have been 
discovered indeed. 

For the sake of completeness let us recall below different criteria for the 
isochronicity of periodical solutions of equation (1). 

Using formula (4), Landau and Lifschitz [9] deduced the following 

Proposition 3-1 [9, Chap. 3] When g(x) is continuous and hypothesis 
(H) holds, (1) has an isochronous center at the origin if and only if 

x-A(x) = 2^2G(x) 

for all < x < b where A is the involution such that G(A(x) = G(x) and 
A(x)x < for x 7^ 0. 

When g(x) is a continuous function, g{x) and x having the same sign, 
Koukles and Piskounov [7] produced necessary and sufficient conditions so 
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that the center of the system (2) is isochronous. 



Proposition 3-2 [7, Th 5] A set of necessary and sufficient conditions 
for the period of every solution of (1) near to be equal to a constant T is: 

1 - g(x) is continuous and positive for small positive x. 

2 - liminf x ^ | ^ M 0. 



3 - T > limsupx^Q '" 



\/g( x )/ x 

4 - g{—x) = — ^{^\/2x — (Jq g^du)^ 1 ]^ 1 where index —1 denotes an in- 
verse function. 



This result was improved by Koukles and Piskounouv themselves in the 
analytic case. They proved the following 



Proposition 3-3 [7, Th 6] Let g{x) be a real analytic function. Then 
the center of the equation (1) x + g(x) = is isochronous if and only if 
the inverse function x — Q(z) 

J 

is of the form 

where P is a real analytic function such that P(0) = . 

Later, Urabe proposed some refinements of Proposition 3-3 by considering 
the assumption of the differentiability of g(x) at and proved the following 
result which is most used than Propositions 3-2 and 3-3 

Proposition 3-4 [13] Let g(x) be a C 1 function defined in V a 
neighborhood of verifying xg(x) > in V /{0}. Then the system (2) has 
an isochronous center at the origin if and only if g(x) may be written 

9{x) = T^hjx) 

where h(X) is a C l odd function and X = ^2G(x), ^ > for i^O. 
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We will propose an alternative approach in order to derive isochronous 
potentials. Some other criteria or equivalent characterizations will be pre- 
sented. Their significance makes the study of the isochronicity much easier, 
since any isochronous potential G appears to be solution of a differential 
equation. More precisely, we state the following 

Theorem B Suppose hypothesis (H) holds and let g{x) be an analytic 
function and G{x) = f£ g(s)ds. Then the equation 



where f is an analytic function defined in some neighborhood ofO. 

We first deduce from Theorem B that the analytic functions / and g are 
naturally related and by (10) one has necessarily 



x + g(x) = 



(1) 



has an isochronous center at if and only if 



d 

dx 



[G{x)/g\x)\ = f(G) 



(10) 



/(0) 



/(0)/3,f (0) = -(7/%" 3 (0) + <? (4) (0)/5, 



f(0) 




1 



( 6 )(0) 



310 
~54~ 



/ 5 (0) + 2/ 2 (0)s (4) (0),.... 



Proof of Theorem B Recall that 




and by Lemma 2-1 its derivative T' 



(c) = 4jL may be written 




(11) 




(12) 
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Suppose at first that G R 2 is an isochronous center of system (2). Then, 
by Proposition 3-1 



x - A(x) = 2^2G(x) for all < x < b. (13) 
Deriving this expression 

1 - A'(x) = 

it implies 

G(x) 2 



g 2 (x) (1-A(x) 2 ' 
Therefore, by deriving with respect to x 

d,G(x)^ d 2 4A"(x) 



dx v g 2 (xy dx\l-A(x) 2i (1 - A'(x)f 
On the other hand, we need the following 



(14) 



Lemma 3-1 For any analytic involution A(x) defined for all x G [a, b] 
the following expression holds 

A"(A(x)) A"(x) 



(1 - A'(A{x))) 3 {1-A'{x)f 
Proof Indeed, to prove this lemma we derive A(A(x)) = x it yields 



■\'(M.r)) -i- and A"((A(x)) A " U) 



A'(x) vv v " A' 3 (x) 

We now replace x by A(x) in the right side of (9), one obtains 

A"(A(x)) A"(x) 1 

(1 - A'(A(xW ~ A«(x) (1 - ^ 

So, since A'(x) ^ one gets after simpification 

4A"(A(x)) 4A"(x) 
(1 - A'(A(x))f ~ (1 -A'(x)f 
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This lemma implies that expression (14) is A-invariant and therefore 
J^[p^y] is an analytic function only dependent on G since G(A(x)) = G(x). 

We have now to prove the converse. We will use for that the following 

Lemma 3-2 The derivative of the period function T(c) may be written 

as 



T'(c) = (— =) / ' d ; V " g{x)dx. 

cV2 Jo ^c-G{x) 



Proof Indeed, we have seen that 

d :(G/g 2 )(x) 



c - G(x) 



g(x)dx 



for a < a < and < b < b. 
By splitting the integral we get 



cV2 A ^Jc-G{x) cV2 ^ c _g(x) 



+ ( — -=) / ^ =g(x)dx 



Thus, by this lemma condition (10) J^[G(a;)/g 2 (:r)] = /(G) implies nec- 
essarily that £[G(x)/g 2 (x)] = £[G(x)/g(A(x)) 2 ]) and then T'(c) = 0. 
Theorem B is then proved. 

Let us consider F the primitive of / such that F(0) =0. Then 

(10) j-[G{x)/g 2 {x)\ = f{G{x))^g{x)j-[G{x)/g\x)\=g{x)f{G{x)). 

Integrate by parts, it yields |g - x = F(G(x)). Therefore (10) is equiv- 
alent to 2G(x) — xg(x) = g(x)F(G(x)). Thus, we get another criteria of 
isochronicity 
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Corollary 3-3 Under hypotheses of Theorem A, equation (1) admits 
an isochronous center at if and only if 

2G{x) - xg(x) = g(x)F(G(x)) (10') 

where F is an analytic function defined in some neighborhood of 0. 

Some other equivalent conditions may also be deduced 

Corollary 3-4 Under hypothesis (H), is an isochronous center of 
(1) if and only if x = x{G) is an analytic solution of the linear ODE 

where f is an analytic function. 

Moreover, this solution must satisfy the conditions: 

x(0) = 0, lim G ^ (—) = I. 



Proof of Corollary 3-4 Let us consider again F(t) the integral of f(t). 
Then, it is easy to see that condition -^(G/g 2 ) = f(G) is equivalent to 

2G^ = x + F(G). (16) 

We derive (16) with respect to the variable G, we then obtain (15). 
So, for any analytic function / the linear equation (15) admits a unique 
solution x = x(G) according to initial conditions. More precisely, consider 
the change x = \/2G + y. Then y = y(G) is solution of 

with initial points : y(0) = 0,y'(0) = 1. A resolution of the last equation 
yields 

rG F (v) 
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Thus, a solution of (15) may be written 



, — r G F(u) 



(2l/) 3 /2 



Let us consider now the involution A such that G(A(x)) = G(x). Another 
consequence of Theorem B is 

Corollary 3-5 Let G(x) = Jq g(s)ds be an analytic potential of the 
scalar equation 

x + g(x)=0. (1) 

Let A(x) be an analytic involution defined by: 

G{A{x)) = G(x) and A(x)x < 0. 

Then is an isochronous center of (1) if and only if 

4-[G(x)/g*(x)]= [A " U) 



dx L ^ "* ^ » (1 - A'(x)f 
Moreover, the last expression is A-invariant. 



Proof of Corollary 3-5 Let A be an analytic involution then 
G(A(x)) = G(x) implies 

f x=g{x) = i A g{Mx ) ) 



and 



-M{x)) = {^f%{x). 
g z ax g z 



Deriving the last expression we then obtain 

A'(x) d ^-(A(x)) = 2A'(x)A"(x)%(x) + A' 2 (x)^-(x). 
Since A'(x) ^ then, after simplification we get the differential equation 

d -^-(A{x)) = 2A"(x)^(x) + A'(x) d ^-(x). (17) 



By Theorem B, 

£ [G{x)/g\x)\ = /(G) implies ±[G{x)/g\x)\ = ±[G(x) / g(A(x)f]). 
Therefore, the solution of equation (17) is 

?<*> = i^hw (18) 

since A'(0) = -1. Thus, by deriving one gets £[G (x) / g 2 (x)] = (1 ^^ )3 . 
To prove the converse we require again Lemma 3-2 

l r b j- (40*0) - f (4(^0*0) 

T'( C ) = (-L) / ^* 21 7 d ; g 1 ^ (x)^. 

cV2 Jo ,Jc-G(x) 

The analytic involution A(x) can also be defined as a solution of a linear 
ODE. The following result is analogous to Corollary 3-4. 

Corollary 3-6 Under hypotheses of Theorem A, equation (1) admits 
an isochronous center at if and only if the involution A = A{G) is a solution 
of 

HA 

2G— = A(G) + F(G) (19) 

/ is an analytic function and F is its integral. 
Moreover, this solution must satisfy the conditions: 

A(0) = 0, Zim G _ M) (^) = 1 



Proof of Corollary 3-6 To see that, we start from A(A(x)) = x and 

dx 
dG 



G(A(x) = G(x). So equation (11) 2G% = x + F(G) is equivalent to 



2G (A(x) = dG{ f x)) [A(x) + F(G(A(x))} = ^-(x + F(G)). 
ax ax 

On the other hand, deriving G(A(x) = G(x) with respect to x we get 
(^(dG^xii) = d^i_ We then obtain 

A{x) + F{G{x)) = [x + F(G(x))] dAU ' 



dx 
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Hence we deduce the ODE (19). 

Remark 3-7 The method of Urabe requires the use an intermediary 
function h that is not in general explicitly known. Indeed, in order to de- 
termine the potential isochrones one must be able to show by means of the 
change of variables X = \[2G that the inverse function of G(x) = y is of the 
form x = X + H(X). What is not easy to achieve in any case. On the other 
hand, our approach is more direct. Our criteria are simply relations between 
an isochronous potential and its derivative. To be clearer, condition (10) 
of Theorem B is a sort of non linear differential equation of order one with 
respect to the variable x. All solution G = G(x) of (10) j-[G(x) / 'g 2 (x)] = 
f(G(x)) or equivalently (10') 2G(x) — xg(x) = g(x)F(G(x)) provides an 
isochronous potential verifying G(0) = G'(0) = 0, G"(0) = 1. 

It is well known that the harmonic potential G = \x 2 is the only polyno- 
mial potential which is isochronous. Concerning the rational case, Chalykh 
and Veselov [2] proved the following 

Proposition 3-8 ([2]) Under hypothesis (H) a rational potential G{x) 
(which is not a polynomial) is isochronous if and only if 

G(x) = ^- 2 [ax + l-^—r. 
8cr ax + 1 

The proof we give below is different of that given by [2]. 

Proof of Proposition 3-8 Let us recall that for a involution A we get 
G(A(x)) = G(x) 

and A(x)x < for all x G [a, b]. By Proposition 3-1 this potential is 
isochronous if 

x - A(x) = 2^2G(x) for all 0<x<b 
or equivalently it verifies the functional equation 

G(x) = G(x - 2^2G{x) for all < x < b. 

When the potential is rational it is required that G(x) has to be the square 
of a rational function of the form 

cm = \m? 
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where P(x) and Q(x) are polynomials without common zeros. 

This means the involution has to be meromorphic and has the same poles 
as G(x). It can be written 

2xP(x) xQ(x) - 2xP(x) 

[x) " x ~ [ ~oW ] " W) 

But, it is known (see for example Theorem (15.4) p. 296 of [6]) that the 
only meromorphic functions verifying A 2 (x) = x are A(x) = L~ l (jj^) where 
L(x) = ax + b a affine function since other meromorphic functions are not 
invertible. A is then an homographic function 

Thus, by hypothesis (H) and since A(0) = and A'(0) = — 1 we get 

A(x) = X — 

implying 

2P(x) 1 ax + 2 
— = 1 H = . 

Q(x) ax + 1 ax+1 



4 Parametrization of isochronous centers 

In this part we suppose that all functions are analytic. In particular, let us 
write 

g ( x ) = x + J2 a nX n and G(x) = -x 2 + J2 -^x n+ \ 

n>2 ^ „> 2 U + 1 

and suppose that r is the radius of convergence of these power series. By 
the Cauchy-Hadamard formula ^ = livrin^^ \ a n \ l l n . 

The purpose of this section is to highlight conditions so that an analytic 
potential G(x) be isochronous. To that end let us write 

G(x) = \x 2 + Gi(x) + G 2 (x) 

where the function G\(x) = J2k>2 ^zf^^ 1 i s °dd and 
G 2 (x) = J2 k >2^x 2k is even. 
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We looking for general conditions on coefficients a n ensuring the isochronic- 
ity of the center of equation (1) x + g(x) = 0. Recall that only the har- 
monic potential G(x) = \x 2 is an even isochronous potential [G\[x) = 0). 
That means the odd coefficients a 2 k+i cannot be all zero when the isochronous 
potential is non-harmonic. In fact, we will prove little more. 

Starting from Theorem A we will show there are infinitely many neces- 
sary conditions verifying by the coefficients in order the potential G(x) be 
isochronous. More precisely, we state that the even coefficients a 2 k may be 
free and the odd coefficients a 2 k+i are polynomials with respect to a 2 k- 

Theorem C Let the analytic potential 

G(x) = lx 2 + J2 —% n = + G^x) + G 2 (x) 

2 n>3 H 2 

of equation 

(I) x + g(x) = 0. 

When the equation (1) has an isochronous center at then the odd coefficients 
of the expansion of g{x) can be expressed in terms of rational polynomials 
involving the even coefficients: 



In particular, when the potential G(x) is isochronous thenG 2 {x) = is equiv- 



alent to G\{x) = 0, i.e. G{x) = \x 2 is harmonic. 



Thanks to Maple we are able to calculate the first terms : 

10 2 14 56 4 

a 3 — -^a 2 i o 5 — — a 2 a 4 — ^ a 25 

-592 3 848 6 24 36 2 

a 7 = -^-a 4 a 2 + —a 2 + —a 2 a 6 + — a 4 , 



110 440 3 27808 5 536800 8 1144 2 2 22 

a 9 = ^a 2 a 8 - ^a 2 a 6 + ^^-«2 a 4 - 65gl a 2 - ^-a 2 «4 + ya 4 a 6 
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„ _ 52 , 57616 n 4 n 2 2600 _ 3 „ i 125008 „ 5„ 4837664 7„ , 

Oil — YT °2 + ^35~ a 2 a 4 sF ° 2 fls 56T - 2 6 3645 a 2 a 4 + 

5631808 _ 10 2392 3 7384 2 _ _ , 52 , 78 _ 2 

6561 02 _ "125" 0(2 ^ ~ I05 a 2 a 4 a 6 + JE a 4 a 8 + a 6 



ai3 — —72 ag a§ a^ 2 — ^p 2 - ag 2 ^ a§ + 38 2 1 7 76 ag 4 ^ + y§ azs + ff a 4 a io ■ 

in/Q„ „ 9430624 n 2 „ 6 1 375769408 _ _ 9 1 166544 „ 3 „ 3 920 „2„2 

iu/o a§ 1215 H ig68 3 °4 a 2 ~\ 225~ a ^ 2 — 21 2 0(5 ' 

2190080 „ 7„ 14000 „ 3„ _i_ 300944 „ 5 _ 74681600 - 12 616 „ 4 
729 ^ fle ~ ~29T 0,2 ai ° + ~W ^ 0,8 656T~ 0,2 ~ 125 a 4 



Proof of Theorem C Let the analytic function 

n 

g {x) = x + akX k - 

k=2 

By Theorem A the potential G(x) = \x 2 + J2n>3 ^ 1 x n is isochronous if and 

only if the following equality holds (10) ^r(p^) — f{G) where / is an 
analytic function, set 

/(G) = b + hG + b 2 G 2 + b 3 G 3 + .... 
After replacing and equaling the two sides of (10) : 

(Ti \ Tl O / Th ft 1 \ / ft \ ^ 

* + -(, 2 +E^ +1 ) 1 + E ^ x + E^M = 

fc=2 / fc=2 K "+" 1 \ fe=2 X / V k=2 ) 

Then we identify the analytic expansions of the two expressions. The 
unknown coefficients will then be determined by comparing powers in x: 

2 10 2 2 24 

»0 = -o°2, °3 = 7T«2 & 1 = lla 2a 3 - 6a 2 - — d 4 

o y 

116 

-(10/3) a 5 - (1/3) 6ia 2 + — a 2 a 4 + (14/3) a 2 4 + 4 a 3 2 - 13 a 3 a 2 2 = 0, 

15 

30 1 1 , 20 = 21 2 35 2 3 

— z- a 6 -- 6 2 -- Dia 3 +10a 2 a 5 — — a 2 +— a4a 3 -17a 4 a 2 — — a 2 a 3 +25 a 3 a 2 =0,.... 

From these recursion formulae the coefficients can be easily determined. 
After eleminating b ,bi,b 2 , ... we then deduce the expressions of coefficients 
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°3? °5) a li ■■■■ 

The last part of Theorem C will proved by recurrence. Suppose that all 
coefficients a n = for any n <2p where p is a positive enteger and bo = b\ = 
b 2 = .... = b p -i = 0. So, we may write 

g(x) = x + a 2p x 2p + a 2p+i x 2p+1 + a 2p+2 x 2p+2 + .... 

and 

G (x) = -x 2 + -^^x 2p+1 + -^±1^+1 + ^E+^ X W + ... 
K ' 2 2p+l 2p + 2 2p+l 

Calculate 

^(^) = -2pa 2p x 2p ~ l - a 2p+1 x 2p (2p+l) + a 2p x 2p+l + a 2p+1 x 2p+2 + 

and equaling with 

f(G) = b p G p + b p+1 G p ^ + ... = K 2p + |±1^+ 2 + |±f^+ 4 + ... 
Thus, necessarily a 2p = 0. 

Remarks 4-1 

Consider again Theorem B then equality 

(10) J-( G A? 2 ) = b ° + b ^ G + ^G 2 + 6 3 G 3 + hG 4 + ... 
ensures the isochronicity of the potential G. 

It is also possible to express coefficients a n of the analytic expansion of g(x) 

(or its successive derivatives at 0) in terms of 6 , b 1: & 3 , .. 

Indeed, 

g'(0) = 1, <f(0) = -3b ,g"'(0) = 15^,^(0) = -1056^ - 30^, 
^(5)(o) = 945^ + 6306 &i,# (6) (°) = -103956^ - 11340^^ - 8406 2 , 
gW(0) = 2079006[]6i + 302A0b b 2 + 1351356^ + 113406?, ... 

g may also be written 

g(x) = x-l b x 2 +l 6oV+(-^ h - f b 3 ) x 4 +{- 5 b (f* b 3 + f h) x 5 + 
(- IS 6 « - f " 1 &o 5 ) * 6 + (f + f C + ^ + £ &i 2 ) a: 7 + ... 
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and the equation x + g(x) =0 has an isochronous center at for param- 
eters values b ,bi,b 2 , ... 

Thus, any isochronous potential G is a multiparameters analytic function 
which may be expressed under the general form 

G(x) = \x 2 -\ b x 3 +l 6oV+(-i h - | b 3 ) x 5 +£- b (fi 6 3 + f 61) x 6 + 

if 7 t 21 l 27 231 l 5\ 7,1 ( 55 u 3 h , 429 L 6 , 3 L L , 1 L 2\ 8 , 
7 I, - 120 a _ ~8~ 1 16" ) + 8 V"8" °° ° l + T6" + 10 °°° 2 + 16 1 / + 

A general expression of a 2p as a polynomial of b , b±, b 2 , ... seems difficult 
to built. Nevertheless, it is possible to know the first coefficient 

-1 2p + l 
flw "2H(2p)(2p-l) Vl + - 

5 Application to the search of isochronous po- 
tentials 

In this part, we will see that above results allow us to determinate families 
of isochronous potentials. More precisely, we then apply Theorems B to pro- 
duce potentials with constant period. Some of them are new 

A three-parameters family of isochronous potentials 

To be concrete, we derive at first a three-parameters family of potentials 
which appears little more general than the one given by Dorignac [6]. 

Let us consider the following case 

d ,G . a 
^V )= (1 + /3G) 3/2 ' 
where a and j3 are real parameters such that 2a 2 < (3. 

Thanks to Maple a resolution of these equations yields 
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G(x) 



8a 2 + + 2a 2 )(4ax + (3x 2 ) - (4a 2 + 2af3x)J2(2 + (5x 2 + Aax) 



2(f3-2a 2 ) 2 

(20) 

Then, the above potential is isochronous according to Theorem B. It may 
also be writen 



G(x) = \x 2 



where 



2a + /3x - a^J2{2 + (5x 2 + Aax) 
X= p -2a 2 

The function h of Proposition 3-4 is 

HX) ^ x 



V2TW 1 ' 

Recall that X is such that ^ > so that X is a bijection of x. Moreover, 
h(X) is odd function verifying | h(X) |< 1 for any X e R since 2a 2 < (5. 

Applying scaling property of isochronous potentials, (Claim 2, Corollary 
of [6]). The potentials G(x) and -^G^x) have the same period. That means 
the following three-parameters potentials family is isochronous 



1 „ [2a + /3 7 x - a x 2 2 + /3 7 2 x 2 + 4a 7 a; 2 

G(x) = —X 2 ( 1X ) = — l-i -4 — - (21) 

K J 2 7 2 w ' 2 7 2 (/3-2a 2 ) 2 V ' 

The derivative of G(x) is then 



2a + frx - «[^2(2 + + 4. 7 *~)P 7 - 

9[X) 7 2 (/3-2a 2 ) 2 

The involution A defined by G(A(x) = G(x), A(x)x < may be written 

i [2a + (3-fx - a v /2(2 + /3 7 2 x 2 + 4a 7 a;)] 



A(x) = x - 2- 



7 (/3 -2a 2 ) 

26 



As special cases we may derive the following 



1 - For any 7 7^ and f3 = 2a we obtains the two-parameters family 
of isochronous potentials (see (24) of [6]) introduced for the first time by 
Stillinger and Stillinger [12]. Indeed, after replacing in (21) (5 by 2a and 
simplifying by 2a it yields 



[1 + -yx — J I + a^ 2 x 2 + 2a^x)} 2 

G <*> = V I)- (22) 

2 - The case a = (5 = and 7 = 1 yields the harmonic potential : 
G(x) = \x 2 . 

3 - The case ft — and 7 = 1 gives the Urabe potential (see Corollary 
2-2) : 

a 2 a a 

4 - The case 2a = (3 and 7 = 1 yields the Bolotin-Mc Kay potential (see 
[6]). Indeed, 

d ,G . ct 
^ V ~ (1 + 2aGf 2 

implies 

[1+x - yl + ax 2 + 2ax)] 2 



G(x) 



2(1 -a) 2 

We then deduce the function h of Proposition 3-4 

h(X) aX 



a 



5 - The case 2a 2 = f3 and 7 = 1 yields the isotonic potential (see [6]) 
Thanks to Maple a resolution of 

d ,G , a 



dx K g 2 ' (l + a 2 G) 3/2 

gives 

1 x 2 {2 + axf 
~ 4 (1 + ax) 2 
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The function h of Proposition 3-4 is 

h(X) = aX 



VI + X 2 a 2 
and its integral is 

H{x)= VTTx^_i_ = x _ x 

a a 

Thus 

We then obtain by another way 

4 (1 + ax) 2 8a 21 ax + V 

Others isochronous potentials 

We will give others two-parameters families of potentials with constant 
period which seems to be new in the literature. 

1 - Let us consider 

d ,G . a 



dx K g 2> (1 + 2/3 2 G) 5/2 ' 
Then the function h of Propositon 3-4 may be calculated 

h(X) = a ( 1/3 ^—3^ + 2/3 ■ 



[1 + X 2 (3 2 f 2 1 VT+xW 

and its integral is 

a{l + 2X 2 p 2 ) a 



H{X) = 1/3 



(3 2 y/TTXW 3/3 2 



Taking for example (3 2 = 2a/ 3 and thanks to Maple we find another 
isochronous potential. 
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Indeed, since x = X + H(X) and solving H(X) 2 — (x — X) 2 we get 

= (1/24) + (1/24) ( " A ' + ^-g ),3+4 *° ) + (1/24) ^ 

where 

u(x) = - 5751 + 31536 xa + 21600 a: 2 a 2 + 768 x 3 a 3 + 256 x 4 a 4 

+72 \/3\/ 1024 :r 6 a 6 + 4608 x 5 a 5 + 26496 a; 4 a 4 + 62208 x 3 a 3 + 101412 x 2 a 2 + 86022 xa + 45927 



So, according to Theorem A and in using scaling properties the potential 

[ttfrs)] 1 / 3 , (l6f 2 x 2 a 2 +24 1 xa- 
v(3+4 72:a) afu^x)] 1 



1 3+47jra]2 

U W — 1 01(3+472:0:) ^ o[«(7a;)l 1 /3 1 o J 



is isochronous. 
2 - Let 



d G 2a + 2a 3 G 



dx y g 2> (1 + 2a 2 G) 5/2 ' 
The function h is then 

. , v . _ 2aX + a 3 X 3 
k[X) ~ (l + a 2 X 2 ) 

and its integral has the following simple form 

H(X)= ^ 



VTTc^x 2 ' 

Thanks to Maple we find 

x _ 1 (-10aV + 1 + a 4 x 4 ) l + a¥ 

Qa 2 x Qa 2 xv(x) 6a 2 x 

where 

v {x) = [-15a V + 39a V + 1 + a 6 x 6 + a 3 x 3 yJ(-3 + 33a 2 x 2 + 3a 4 x 4 )] 1/3 
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